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In stellar core-collapse events matter is heated and compressed to densities above nuclear matter 
saturation density. For progenitors stars with masses above about 25 solar masses, which eventually 
form a black hole, the temperatures and densities reached during the collapse are so high that a 
traditional description in terms of electrons, nuclei, and nucleons is no longer adequate. We present 
here an improved equation of state which contains in addition pions and hyperons. They become 
abundant in the high temperature and density regime. We study the different constraints on such 
an equation of state, coming from both hyperonic data and observations of neutron star properties. 
In order to test the zero-temperature versions, we perform numerical simulations of the collapse of a 
neutron star with such additional particles to a black hole. We discuss the influence of the additional 
particles on the thermodynamic properties within the hot versions of the equation of state and we 
show that in regimes relevant to core-collapse and black hole formation, the effects of pions and 
hyperons on pressure, internal energy and sound speed are not negligible. 

PACS numbers: 26.60.Kp,26.50.+x 



I. INTRODUCTION 

Supernovae and hypernovae figure among the most 
spectacular events observed in the universe because of 
the immense amount of energy involved. In general, one 
can distinguish between thermonuclear and core-collapse 
events. Here we shall be interested in the latter. These 
occur at the end of the life of massive (M > 8M0) 
stars: when the iron core exceeds the Chandrasekhar 
mass a gravitational collapse is induced. At the cen- 
ter a compact star is formed, which is a neutron star in 
the classical gravitational supernova. Depending among 
others on the progenitor mass, its metalicity, and rota- 
tion, as well a black hole can be formed. These events 
are known as hypernovae or collapsars. For about thirty 
years, simulations have been performed in order to ex- 
plore these events and to answer related questions, for 
example on the precise conditions for forming a neutron 
star or a black hole. The simulations are extremely com- 
plex, since they involve many different ingredients: multi- 
dimensional hydrodynamics, neutrino transport, general 
relativity and complicated microphysics. Despite all the 
effort, many unknowns remain in the simulations, in par- 
ticular on the engine driving a successful supernova ex- 
plosion. Apart from the observations via electromag- 
netic radiation, the neutrino and gravitational wave sig- 
nal could give interesting information on the models. 

The microphysics input for the simulations concerns 
essentially two domains, the rates for neutrino-matter in- 
teraction and deleptonization, i.e. electron capture, and 
the equation of state (EOS). Here we will discuss the 
latter. It is not an obvious task to construct an EOS. 
The main difficulty arises from the fact that very large 
ranges of (baryon number) densities (10~^° fm~^ < hb ^ 
1 fm~'^), temperatures (0 < T < 150 MeV) and hadronic 



charge fractions [Q <Yq — Uq/ns ^ 0.7) have to be cov- 
ered. Uq here denotes the total hadronic charge density, 
which in many cases is just given by the proton density. 
Within this range the characteristics of nuclear matter 
change dramatically, from an ideal gas of different nu- 
clei up to uniform strongly interacting matter, contain- 
ing in the simplest case just free nucleons but potentially 
other components such as hyperons, nuclear resonances 
or mesons. Even a transition to deconfined quark matter 
cannot be excluded. Although there is a large variety 
of EOSs available for cold dense matter relevant for the 
description of neutron stars (see for example [1] and refer- 
ences therein), at present, only a few hadronic EOSs ex- 
ist which are commonly used in core collapse simulations, 
where temperature effects play a crucial role. There is the 
one by Hillebrandt and Wolff [2], used by some groups 
performing supernova simulations, that by Lattimer and 
Swesty [3] as well as that by H. Shen et al. [4]. The 
two latter, publicly available, are most commonly used 
in core-collapse simulations. They use different nuclear 
interactions, but are based on the same limiting assump- 
tions: they take into account non-interacting a-particles, 
a single heavy nucleus and free nucleons in addition to 
the electron, positron and photon gas. 

However, in particular at low densities, i.e. below 
roughly nuclear matter saturation density ns ng 
0.16 fm^"^ (corresponding to a mass density of about 
^ 10^'* g/cm^)^ the composition of matter is much more 



^ We shall work here exclusively with baryon number densi- 
ties, since the baryon number is a conserved quantity, notably 
throughout a hydrodynamic simulation, contrary to the mass 
density which is not conserved. Many codes, for dimensional 
reasons, work, however, with a mass density. The latter can be 
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complicated, with a large number of different nuclei. Al- 
though this shoiild not have a large impact on the purely 
thermodynamical properties [6], it is important to cor- 
rectly describe the composition of matter in order to de- 
termine the electron capture rates and neutrino interac- 
tions. Therefore, in the last years, several groups have 
started to build EOSs using mainly statistical approaches 
to improve the low density part of the EOS (see e.g. [7 
13]). It has been shown that especially the presence of 
additional light nuclei can have an influence on the su- 
pernova dynamics and among others on the neutrino sig- 
nal [5, 11, 14, 15]. We do not discuss this point in the 
present paper since we are mainly interested in a high 
density and high temperature extension of the EOS. 

Up to know, there are less attempts to improve the 
high density (ns ^ no) and high temperature (T ^ 20 
MeV) part of the EOS, although there are many indi- 
cations that probably the physics of the standard EOS 
is too poor in this regime, too. First of all, our knowl- 
edge about the quantum chromodynamics (QCD) phase 
diagram suggests a transition to the quark-gluon plasma 
(QGP) within the range of densities and temperatures 
reachable in core collapse events, i.e. within the range 
of our tables. Of course, there are lots of uncertainties 
about this phase transition, so that its occurrence cannot 
be affirmed, but the possibility has to be kept in mind 
when employing a purely nuclear EOS such as the two 
EOSs by Lattimer and Swesty [3] (LS EOS) or by Shen 
et al. [4] (Shen EOS) up to densities well above nuclear 
matter saturation density and temperatures as high as 
several tens of MeV. There is some first work including 
this phase transition, see [16]. Secondly, even without 
thinking about a QCD phase transition, other forms of 
(non-nucleonic) matter should appear at high densities 
and temperatures. Already for a long time for cold EOS 
used for neutron star models, hyperons, pions and kaons 
have been considered. At temperatures above about 20 
MeV, this point becomes even more crucial. This has 
been confirmed by the first attempts to include hyperons 
and pions in the Shen EOS [4] for simulations, see [17 19]. 
The effect of these high density and high-temperature ex- 
tensions of the EOS on the simulations is not negligible, 
sec e.g. [16, 17, 20, 21]. Let us in particular mention that 
Sagert et al. [16] found that the QCD phase transition 
could induce a second shock wave which in their simu- 
lations leads to a successful explosion. We shall discuss 
here the construction of a new EOS including hyperons 
and pions based on the Lattimer and Swestj^ [3] (LS) 
EOS and the effects on some thermodynamic quantities 
important for the simulations. 

The paper is organized as follows. In Section II we 
briefly recall the basics of the LS EOS [3] upon which our 



obtained easily just by multiplying the number density by a con- 
stant mass, e.g. the neutron mass ran, see the comment on this 
point in [5], too. 



model is based. In Section III we present our extension 
including hyperons and pions. In the following Section, 
Section IV we discuss the existing constraints we have on 
the construction of the extended EOS. In particular we 
discuss the compatibility with the recent observation of a 
neutron star with almost two solar masses [22], claimed 
to exclude the existence of additional particles such as 
hyperons, mesons or quarks within cold neutron stars. 
Section V gives an illustration of the usability of the hy- 
peronic EOSs at zero temperature and beta equilibrium. 
Sec. VI is devoted to a discussion of the results at finite 
temperature and we conclude in Section VII. 



II. THE LATTIMER AND SWESTY EQUATION 
OF STATE 

Let us start the description of our model for the ex- 
tended EOS with a description of the original EOS by [3] . 
We have chosen this EOS as basis for our work in order to 
have an approach for the hadronic interaction at hand dif- 
ferent from the attempts to include hyperons in the Shen 
EOS, employing a relativistic mean field model [17-19]. 
The motivation is of course that there are large uncertain- 
ties on the hadronic interaction, so that it is interesting to 
compare two different types of models. In addition, the 
LS EOS is one of the two commonly used EOSs in com- 
putational astrophysics, so that a comparison of existing 
results in the literature with results from our extended 
model should be simplified. 

As mentioned above, the LS EOS [3] models the mat- 
ter as a mixture of one (average) heavy nucleus, a par- 
ticles, free nucleons, electrons, positrons and photons. 
Electrons and positrons are treated as non-interacting 
relativistic gas in pair equilibrium, neglecting electron- 
screening effects; photons are treated as an ideal ultra- 
relativistic gas. Equilibrium with respect to strong and 
electromagnetic interactions is supposed, while no ^ equi- 
librium is assumed, as expected during core-collapse su- 
pernova. 

Concerning the nuclear part, the LS EOS follows the 
works by Lattimer et al. [23] and by Lattimer and Raven- 
hall [24]. Some simplifications have been made with 
respect to Refs. [23, 24], e.g., the neutron skin is ne- 
glected and a simpler momentum-independent nucleon- 
nucleon interaction is employed instead of a standard 
non-relativistic Skyrme parameterization. Within the in- 
homogeneous phase at low density, nuclei are supposed 
to arrange themselves in a body centered cubic lattice 
which maximizes the separation of ions. According to the 
Wigner-Seitz approximation, each ion is at the center of a 
neutral-charged cell, surrounded by a gas of free nucleons, 
a and electrons. Interactions between the outside gas and 
the nuclei are taken into account through an excluded 
volume. Nucleons are treated as non-relativistic parti- 
cles; a-particles as hard spheres of volume Va = 24 fm^ 
forming an ideal Boltzmann gas. As the density in- 
creases, nuclei undergo geometrical shape deformations. 
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until they dissolve in favor of homogeneous nuclear mat- 
ter above approximately saturation density. The forma- 
tion of non-spherical nuclei (^^pasta-phase") is described 
by modifying the Coulomb and surface energies of nuclei, 
as discussed in Section 2.8 of Ref. [3]. The phase tran- 
sition to hulk nuclear matter is treated by a Maxwell 
construction between the two phases. The configuration 
of matter and the balance between the different phases 
is given by the most thermodynamically favorable state, 
i.e. the one which minimizes the Helmholtz free energy 
of the system. This procedure, minimizing the free en- 
ergy, guarantees that the LS EOS is thermodynamically 
consistent. 

Let us stress, however, one point concerning the de- 
scription of the transition between homogeneous and in- 
homogeneous matter in the LS EOS. As discussed in 
Ref. [12], it is not satisfactory, since for all subsaturation 
densities matter can be viewed as a mixture of nuclei and 
free nucleons with consequences on the thermodynamic 
properties. In particular, all thermodynamic quantities 
are perfectly continuous. As mentioned already in the 
context of the distribution of nuclei in the inhomogeneous 
phase, our main interest here is a discussion of the high 
density and high temperature part, taking a consistent 
and commonly used EOS for the remaining part. 



A. Characteristics of the Lattimer and Swesty 
equation of state 

The nuclear interaction in the LS EOS contains sev- 
eral parameters, which have been chosen to reproduce 
reasonable values for properties of symmetric (i.e. equal 
number of protons and neutrons) bulk nuclear matter at 
saturation density, for details see the original work [3]. 
These quantities are related to a power-series expansion 
of the energy per baryon around saturation density at 
zero temperature and for symmetric matter: 
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where x = Ub /no — 1 is the deviation of the baryon num- 
ber density from saturation and 13 = {rin — n^jnB = 
\ — 2Yp describes the asymmetry. The properties of the 
EOS are thereby given by the values of the coefficients, 
no, B, K, K', J, L. Of course, this can only give an indi- 
cation on the general behavior of the EOS, since they are 
defined at saturation density and for symmetric matter, 
whereas in the context of neutron stars and core collapse 
events very asymmetric matter at very different densities 
is encountered. 

Nuclear experiments give constraints on the proper- 
ties of the saturation density, no, the binding energy, 
B, the incompressibility, K and the symmetry energy 
at saturation, J. Typical values for ng lie in the range 
0.15 fm^"^ < no < 0.17 fm^"^ and the binding energy is 



15.6 MeV < B < 16.2 MeV. The value of K, roughly 
speaking, determines the stiffness of the EOS, the higher 
the value of K , the stiffer the EOS. But, as mentioned 
above, it is determined at saturation density and for sym- 
metric matter, such that this interpretation has to be 
regarded with caution. Nuclear physics experiments on 
the breathing modes like the isoscalar giant monopole 
resonance give a value for K at saturation density of 
240 ± 10 MeV [25]. The obvious error is rather smafi, 
the result is, however, not uncontested. In particular, 
the extraction of this value from data on isoscalar gi- 
ant monopole resonances depends on the density depen- 
dence of the nuclear symmetry energy, a quantity bin- 
der intensive debate in recent years. We thus think 
that a larger range of values has to be considered. The 
commonly assumed range for the symmetry energy is 
28 MeV < J < 34 MeV. For the other parameters, 
the skewness coefficient K' , the symmetry energy slope 
coefficient L and the symmetry incompressibility Ksym 
data are not really constraining the value. 

The original LS EOS [3] uses no = 0.155 fm-^ S = 
16.0 MeV and J = 28.6 MeV, values in reasonable 
agreement with the constraints^. With the original rou- 
tines, see [26], three sets of boundary and Maxwell con- 
struction tables are provided corresponding to three dif- 
ferent values of the nuclear incompressibility modulus, 
K = 180, 220, 375 MeV. Following the above discussion, 
the two extreme values for K used in the LS EOS are 
in principle disfavored and the preferred parameter set 
for simulations should be that with K = 220 MeV. We 
shall, however, keep the two other sets for two reasons. 
The first one is purely historical: in many simulations 
the parameter set with K = 180 MeV has been used, so 
that for comparison with the existing literature it is in- 
teresting to have this value at hand. The second one is, 
as discussed above, that the narrow range, K = 240 ± 10 
MeV is not uncontested. In that sense, the range of val- 
ues of the LS EOS represents an extreme variation of 
the nuclear parameter sets, i.e. it can give an indication 
about the uncertainties in the simulations, coming from 
the uncertainties on the nuclear part of the EOS and we 
shall in principle keep all three values. 

Finally, for simplicdty we assume that the nucleon effec- 
tive mass is equal to the bare mass: m* = m. In Lattimer 
& Swesty [3], the (density dependent) effective mass term 
is kept in the equations, so that this assumption can be 
generalized. Indeed, mean field theories predict an aver- 
age effective mass m* /m around 0.6-0.8 (see e.g. Bender 
et al. [27] for a review, and references therein). It has 
been shown that the inclusion of a temperature depen- 
dent nucleon effective mass in nuclei, coming from dy- 
namical correlations beyond mean field, may affect the 
core-collapse dynamics [28, 29]. However, as mentioned 



^ The value of J slightly differs from the one given in Ref. [3] 
(J = 29.3 MeV), see Table 2 in Ref. [5] too. 
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earlier, the details of the nuclear part are not the aim of 
the present work, and we shall keep the original version 
of Lattimer and Swesty [3] . 

Note, however, one minor correction with respect to 
the original code. It has been recognized [30, 31] that 
the original LS EOS underestimated the fraction of a- 
particles. The reason is that the a-particlc binding en- 
ergy Ba has to be measured with respect to the neutron 
mass, as all other energies. 

Let us now show the results of this correction to the 

LS routine. In Fig. 1 arc displayed the abundances, 
the entropy and pressure as functions of density, for 
K = 180 MeV, Yq = Yp = = 0.3 and for differ- 
ent temperatures (T = 1, 2, 3 McV). Note that in this 
case the hadronic charge fraction is given by the pro- 
ton fraction and that it is equal to the electron fraction, 
Ye = (rie- —ne+)/nB due to charge neutrality. These are 
typical conditions that can be found in a core-collapse su- 
pernova simulation. We observe, as expected, that in the 
original LS routine the abundance of a particles is un- 
derestimated; as a consequence, nuclei and free nucleon 
abundances are higher at a given density. The pressure 
of the system (upper- left panel) is not very much affected 
by the corrections to the LS routine, since in this density 
range the contribution from Icptons is the dominant one. 
Fig. 1 can be directly compared with Fig. 4 of Ref. [30], 
where the authors plot the results obtained by the origi- 
nal LS EOS and by their 4-species (neutrons, protons, a 
and ^^Mn as representative heavy nucleus) EOS derived 
assuming nuclear statistical equilibrium (NSE). The re- 
sults obtained by our tables and by the NSE EOS intro- 
duced in Ref. [30] agree. Nevertheless, differences have 
to be noticed, especially for T = 1 McV; this could be 
explained by the fact that in their 4-species EOS, the 
authors assume ^''Mn to be the representative heavy nu- 
cleus, while in LS EOS the mean nucleus varies as a func- 
tion of density in order to satisfy the energy minimization 
condition. This affects the relative abundances and the 
macroscopic properties of the system. 

The lowest values of the density in the original rou- 
tines for the LS EOS is Ub = W'^ fm"^. The physical 
reason is that, in principle at low densities and tempera- 
tures below roughly 0.5 MeV, an EOS depending only on 
temperature, baryon number density and charge fraction 
to describe matter in thermodynamic equilibrium is not 
sufficient and a nuclear reaction network has to be used. 
For many purposes, however, a detailed description of 
matter in this regime is not necessary, and it is thus inter- 
esting to have an EOS at hand for this regime. Recently, 
O'Connor and Ott [32] have generated an EOS table in 
which they employ the LS EOS for densities above the 
limiting value of the routines, and, for lower densities, the 
Timmcs EoS [33], under the assumption that matter is 
composed of an ideal gas of electrons, photons, neutrons, 
protons, a-particles and heavy nuclei with the average A 
and Z given by the LS EOS at the transition. We shall 
follow a slightly different approach, see Section III A. 



III. MODEL FOR THE EXTENDED EQUATION 
OF STATE 

A. Low density regime 

At the densities below the limit of the LS EOS, i.e. 

Ub = 10~^ fm~^, at low temperatures matter is com- 
posed of a gas of nuclei and electrons. At temperatures 
above roughly 1 MeV, nuclei are dissolved in favor of free 
nucleons, sec Fig. 2, where the boimdaries between homo- 
geneous and inliomogcneous nuclear matter are shown. 

In this regime, the densities are so low, that the dif- 
ferent particles are only very weakly interacting and a 
description in terms of an ideal gas is completely suffi- 
cient. This is the reason why the choice of O'Connor and 
Ott [32] to employ the Timmes EOS in this regime (which 
is nothing else than an ideal gas of different species) is 
well justified. We follow the same idea, the only differ- 
ence is that we make another choice for the matter com- 
position, in particular for the nuclei present and for the 
matching with the LS EOS at higher densities. O'Connor 
and Ott take one average nucleus with A and Z obtained 
from the average heavy nucleus of the LS EOS at the 
transition density. In this way, it is however not possible 
to describe a potential variation in A and Z of the mean 
nucleus throughout the density range covered by this low 
density extension of the EOS. We therefore use a NSE ap- 
proach, allowing in principle for many different nuclei to 
appear. The setup strongly follows the NSE statistical 
model of Hempel and Schaffner-Bielich [9] , with the only 
exception that we do not consider any excluded volume 
correction because it is not relevant for such low densi- 
ties. In particular, we take into account the Coulomb 
and temperature corrc;c;tions to the binding energies of 
the nuclei in order to obtain reasonable transition tem- 
peratures to homogeneous matter. 

Concerning the matching to the LS EOS at higher 
densities, we have chosen a matching density of Hb = 
5 X 10"^ fm"^. Although the LS EOS reproduces well 
the character of matter at this density, i.e. an almost 
ideal gas of nuclei, photons and electrons, the change of 
the EOS induces small discontinuities in the thermody- 
namic quantities, due to the different treatment of the 
nuclear part. Since the pressure in this regime is dom- 
inated in any case by the electron pressure, the discon- 
tinuity in the pressure is completely negligible. This is 
not the case for the energy density and this is the reason 
why O'Connor and Ott [32] introduce a constant shift 
in the energy density. We judge that the discontinuity 
is small enough, such that this shift, problematic in a 
general relativistic framework, is not necessary. 

B. High density regime 

We have added to the LS EOS pions, muons and hyper- 

ons. For the first two, no interaction has been assumed 
and they have just been added as a free gas, satisfying the 
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FIG. 1: (Color online) Pressure, entropy and abundances (nuclei, a particles, free neutrons and protons) as a function of 
density, for K = 180 MeV, = 0.3 and T = 1 MeV (solid lines), 2 MeV (dotted lines) and 3 MeV (dashed lines). Thick 
black lines correspond to the results of the original LS routine, while thin red lines correspond to the values obtained with the 
modified routine. The differences arise from the correction to the binding energy of a particles. We have used here the mass 
density, defined as pB = m^nB {rUu being the atomic mass unit), for better comparison with the results of Ref [30]. 




FIG. 2: (Color online) On the left panels, boundaries between homogeneous and inhomogeneous nuclear matter (temperature 
versus baryon number density) for different Yq = 0.05, 0.1, 0.2 and 0.4 are displayed. On the right panels, mass fractions 
(free neutrons X„, free protons Xp, alpha particles Xa, and nuclei Xh) as a function of baryon number density are shown for 
different temperatures T = 0.05, 1, 5, 10 MeV, and fixed Yp = 0.2. 
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overall constraint on charge neutrality. There are many 
works considering these additional particles in cold neu- 
tron star cores, see e.g. [34-40]. As already mentioned 
earlier, there is less work in the context of hot and dense 
matter in core collapse events, although the possibility of 
a delayed collapse to a black hole induced by a transition 
to hyperonic matter has been considered in [41, 42] or 
pions (and to less extent kaons, too) have already twenty 
years ago been considered as possible candidates for the 
hot and dense matter in supernova cores, see e.g. [43]. 
The authors of Ref. [43] argue that the presence of pions 
could increase the temperature of the supernova core and 
increase the number of electron neutrinos and thus lead 
to a higher neutrino luminosity in favor of a successful ex- 
plosion. This idea has, however, not been further pursued 
since there are uncertainties about the employed pion- 
nucleon interaction. It is now commonly assumed, that 
there is an s-wave ttN repulsion, reducing strongly the 
number of pions eventually present in supernova cores, 
thus decreasing the effect described by [43]. This repul- 
sive interaction prevents probably pions from condens- 
ing in cold dense matter as discussed for neutron stars, 
too. It is, however, known that the pion gas is one of 
the main components for matter in heavy ion collisions. 
The difference is of course, that for heavy ion collisions 
much higher temperatures (of the order of 100 MeV or 
more) are reached and the baryon number densities are 
much lower than in the core of neutron stars. Core; col- 
lapse events with massive progenitors are situated some- 
what in between, temperatures can reach the order of 
100 McV, but densities of several times nuclear matter 
saturation density are encountered. Here, we mainly dis- 
cuss the effect of pions on the EOS at high temperature 
where the pion gas should be a reasonable approximation. 
Obviously, without interaction, we cannot prevent a 7r~ 
condensate to form below some critical temperatiu-c, de- 
pending on the density, but we consider introducing a 
realistic pion-nucleon interaction for cold dense matter 
to be beyond the purpose of the present paper. 

Concerning the muons we have to mention the fol- 
lowing point. Since lepton flavor conversion via neu- 
trino oscillations is most probably negligible for core col- 
lapse during the first few seconds [44, 45], muon lepton 
number is conserved independently of electron lepton 
number. Therefore in principle a muon fraction, = 
(n^- —n^+)/nB should be added as variable to the EOS, 
in addition to the electron fraction = {n^- — ) / hb, 
temperature and baryon number density. This would, 
however, mean that the simulation codes should evolve 
muon number, too. For the moment this has not been 
done for several reasons. One is of course the techni- 
cal difficulty in adding another evolved quantity together 
with an additional dimension for the EOS. The second 
is that one expects the number of muons to be much 
smaller than the number of electrons due to the mass 
which is higher by about a factor 200, such that the in- 
fluence of the muons should not be very important. With 
these remarks of caution we show in Sec. VI, under the 



assumption that muons are in thermal equilibrium, some 
examples where the densities and temperatures are high 
enough so that the effects of muons are not negligible. 

The contribution of pions, muons, electrons and pho- 
tons to the pressure are obtained from the expression for 
an ideal gas. 
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(2) 



where /3 = 1/T denotes the inverse temperature, /i^ is 
the chemical potential for particle i and E = y^mf + 
is the single particle energy. The upper sign corresponds 
to bosons and the lower to fermions. 

Hyperons are added by extending the model by Bal- 
berg and Gal [46] to finite temperature. This model is a 
nonrelativistic potential model with the contribution of 
the interaction between particles of type i and j to the 
energy density given by: 

£pot{ni,nj) = (1 ^) {aijUtUj + bijUtjU^rij 
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rii + Uj 



+ nf +'n,)) (3) 



where denotes the baryon number density of species i 
and the factor 1/2 has been introduced in order to avoid 
double counting for the interaction between particles of 
the same type, ti represents the third component of the 
isospin of the respective particle, a^ , bij, Cij, and 7ij are 
the parameters defining the interaction. 

The total baryonic energy density is given by the sum 
of the potential energy, the mass energy, 



E 



n,:m,: 



and the kinetic energy. 



2m,: 



(4) 



(5) 



Note that, in order to remain consistent with the LS EOS, 
we do not take an effective baryon mass into account. 

In order to extent the model to finite temperature we 
did not change the structure of the interaction, but we 
only replaced the expression for calculating the densities 
by its finite-temperature version, 



Hi = 



1 



(27r)3 exp [p{E, - /i^)] + 1 



(6) 



Ei = p'^/{2mi) + Ui + nii thereby denotes the single- 
particle energy for particle i. The kinetic energy densities 
are written in an analogous way, 



(27r)3 eMl3{-Ei - fa)] + 1 ' 



(7) 



see Eq. (8) in Ref. [46], too. This simple approach of 
course neglects possible effects of the temperature on the 
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(phenomenological) interaction. Investigating these ef- 
fects is, however, beyond the scope of the present paper. 

The single particle potentials Ui arc obtained from the 
energy density as dspot/drii. The chemical potentials 
for the different particles are obtained from the following 
relation 

Hi = BiHB + QilJiq + L^Hu + L^IJ'i^ + Sii^s , (8) 

with Bi,Qi, Ll^'^ , Si denoting the baryon number, 
charge, lepton number and strangeness of particle i. Note 
that we use the relativistic definition of the chemical po- 
tentials, i.e. the particle rest mass is included in the 
chemical potential. This is the reason why we have added 
the rest mass to the energy Ei, too. We assume that the 
reaction processes involving hyperons, e.g. the dominant 
A production reaction, n + n ^ n + K + A are in equilib- 
rium. In addition, we assume equilibrium for strangeness 
changing (weak) interactions, such that we can take the 
strangeness chemical potential /x^ = 0. 

We have slightly modified the values of parameters for 
the hyperon-nuclcon (YN) and hypcron-hypcron (YY) 
interaction with respect to the work by Balberg and 
Gal [46] in order to be compatible with current experi- 
mental data on hypernuclei, see Section IV A. This model 
has the great advantage that hyperons are added on a nu- 
clear interaction which is exactly the same as in the origi- 
nal LS EOS by [3] , so that an "artificial phase transition" , 
induced only by matching one nuclear model to another 
and which is thus completely unphysical, is avoided. 

IV. CHOICE OF THE PARAMETERS 

This section will be devoted to a discussion of the ex- 
isting constraints on the choice of the parameters for the 
extended EOS. We start with hyperonic data and then 
have a look on cold neutron stars. 

A. Hyperonic data 

In contrast to the nuclear data, hyperonic data are ex- 
tremely scarce, so that there are large uncertainties on 
the hyperonic interactions. Starting with the description 
of the fundamental hyperon-nucleon (YN) interaction, 
it is a long way off having reached the same precision as 
the nucleon-nucleon (NN) interaction, mainly because of 
the very limited amount of scattering data. On the theo- 
retical side, chiral effective field theory calculations have 
been performed, too [47, 48], improving on the reliabil- 
ity of the YN potentials, but still being far from giving 
conclusive results. First results from lattice QCD simu- 
lations of the YN interaction have become available [49] , 
too. 

Concerning the properties of hyperons in dense nuclear 
matter, there are on the one hand many-body calcula- 
tions, starting from the fundamental interaction. In ad- 
dition to the traditional G-matrix calculations [50-52], 



recently a Hartrec-Fock calculation based on a Viowk po- 
tential has been presented in Ref. [53] . The results are in 
reasonable agreement between the different approaches, 
but large uncertainties remain due to the not very well 
known fundamental interaction, see e.g. the discussion 
in Ref. [53]. It has been shown, too, that the inclusion 
of hyperonic three-body forces (YNN) do not strongly 
change the results [54]. We should, however, mention 
that there are very large uncertainties on the hyperonic 
three-body force. The authors of Ref. [54] assume an ef- 
fective phenomenological form similar to the interaction 
in the model by Balberg and Gal [46] we are employing 
here. They, however, limit the strength of the three-body 
force arguing that it should be less important than the 
nuclear one. Here we follow a slightly different philoso- 
phy: we think that the hyperonic three-body force is not 
known well enough to put any constraint on it ad hoc. 
We shall thus limit it only by data and by neutron star 
observations, see next section. 

On the other hand, single-particle potentials Vyn in 
symmetric nuclear matter have been extracted from data 
on hypernuclei. In that way the empirical value for the 
AiV-potential at saturation density (ns = no), V\n ~ 
—30 MeV, has been obtained. This value is in agreement 
with an analysis of (tt", X+) inclusive spectra on differ- 
ent target nuclei [55] and reproduced by most many-body 
calculations, so that it is commonly accepted. Balberg 
and Gal [46] have adjusted their parameters to this value, 
too. For S~, the situation is somewhat ambiguous. The 
observation of a |;He bound state [56] requires an at- 
tractive potential, whereas the analysis of the {tt^,K^) 
inclusive spectra [55, 57] indicate a repulsive potential, 
possibly up to a value of V^n = 100 MeV at satura- 
tion density. Theoretical many-body calculations show 
a large variety of results, too, see e.g. [53], ranging from 
slightly less attractive values as the AN case to strongly 
repulsive values of up to several tens of MeV. Balberg 
and Gal [46] adopt two different versions, one with an 
attractive potential of the same form as for the A, and 
another local potential form giving rise to a repulsive 
potential. We here take the form given in Eq. (3), but 
choosing the parameters in order to obtain a repulsive 
single-particle SiV potential. Concerning the S7V single- 
particle potential, less data are available. Only a few 
events of S-hypernuclei have been observed, so that it 
is much more difficult to reliably fix the well depth of 
the single-particle potential. Balberg and Gal [46] take a 
range Ven = (^20)-(— 25) MeV, whereas newer data in- 
dicate a less attractive potential of Vsn « — 14 MeV [58]. 

Concerning the hyperon-hyperon YY interaction, the 
situation is rather difficult. Early experiments inter- 
preted in terms of production of several double A- 
hypcrnuclei indicate a rather strong attractive poten- 
tial of the order of Vaa ~ — 40MeV [46]. More recent 
measurements [59] are in favor of much lower values, 
Vaa » — 10 MeV. For other hyperons, no data are avail- 
able. We shall take, in view of the faint knowledge, either 
no yF-interaction at all or a universal YY interaction 
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with different values for tfie well depth. For the isospin 
dependent terms we follow Balberg and Gal [46]. 



B. Neutron stars 

In the center of neutron stars densities of several times 
nuclear matter saturation density are reached, so that 
they present an important test for the EOS of matter 
above uq. In contrast to the hot core-collapse environ- 
ment, neutron stars older than several minutes can be 
regarded as cold from the EOS point of view since the 
temperature reached is well below 1 MeV. In addition, 
/3-equilibrium is achieved and neutrinos can freely leave 
the system, so that the EOS in this case is only a function 
of baryon number density. 

In particular observed masses put constraints on the 
EOS. There are a number of precise mass measurements 
from neutron stars in binary systems, for a compilation 
see e.g. Ref. [1]. For a given EOS and a given central 
density, the mass and radius of a non-rotating neutron 
star can be obtained by solving the equations for hydro- 
statical equilibrium together with Einstein's equations. 
In this case of a spherical star it just gives the TOV- 
system. It is known for a long time, that theoretical 
many-body EOS with hyperons predict maximum masses 
of the order 1.4 Mq or below, incompatible with many 
precisely known neutron star masses [60-62]. The re- 
cent precise measurement of PSR 1614-2230 with a mass 
of 1.97 ± 0.04 Mq [22] thus completely excludes those 
EOSs. The authors of Ref. [22] claim that their data ex- 
cludes any type of EOS with "exotic" contributions, i.e. 
other particles than nucleons and electrons. The argu- 
ment is obvious and well known: adding new degrees of 
freedom to the EOS softens it and thus the maximum 
mass decreases. This simple argument is, however, only 
true without interaction. As has already been shown for 
EOS with a transition to quark matter [63, 64], a re- 
pulsive interaction can cure the problem and allow for 
neutron star maximum masses of 2 Mq or even above. 
For hyperons, this seems more difficult, and within the 
microscopic approaches, the origin of the necessary re- 
pulsion at high densities has not yet been found. Re- 
cently different RMF models have been presented which 
successfully reconcile hyperonic matter with PSR 1614- 
2230, see [34, 35, 65-67]. 

Here, we take a phenomenological point of view and we 
choose parameters for the model by Balberg and Gal [46] 
compatible on the one hand with hyperonic data and on 
the other hand with PSR 1614-2230, without looking for 
a deeper understanding of the repulsion. Thus, we shall 
show that it is possible to reconcile the data, in particular 
a neutron star mass of 2Mq with the existence of hyper- 
ons in dense matter, but without answering the question 
why theoretical many-body calculations, predicting the 
existence of hyperons at densities of about 2-3 no, cannot 
reproduce the 2Mq neutron star. The values of the dif- 
ferent single-particle potentials and the maximum mass 




FIG. 3: (Color online) Gravitational mass of non-rotating 
spherical neutron stars as a function of central baryon number 
density for the different equations of state discussed in the 
text. 



of a spherical neutron star for different parameter sets 
are listed in Table I. The first three digits in the name 
of the parameter set thereby indicate the value of the in- 
compressibility modulus of nuclear matter. In Table II 
we list the corresponding parameters for the hyperonic 
interaction. The neutron star masses as a function of 
central density are displayed in Fig. 3 together with the 
constraint from PSR 1614-2230 and the curves obtained 
with the LS EOS for K = 180 and K = 220 MeV. Note 
that pious and muons have only been included for 220pm. 

Parameter set 180BG takes the softest version of the 
original Balberg and Gal paper [46] with two excep- 
tions: we modified the E parameters in order to obtain 
a slightly repulsive V^n and we neglected any YY inter- 
action. This parameter set of course gives a far too low 
neutron star maximum mass, we include it, however, for 
comparison. It can be seen as an extreme case, represent- 
ing a very (too) soft EOS. Set 220BG takes the version 
with the strongest high density repulsion for YN and YY 
parameters from [46] , with, again, one exception: Vsn is 
chosen repulsive. The other three sets show three ex- 
amples of parameters compatible with PSR 1614-2230, 
still giving reasonable values for the well depths at satu- 
ration density. Let us, however, remark that we did not 
find any viable parameter set with an attractive potential 
for Vj:n- For 220g2.8 and 220g3, Vsn is even strongly 
repulsive, at the limit of what is compatible with the es- 
timates discussed in Sec. IV A. For 220pm it is weaker, 
but in that case, no hyperons are present in cold neutron 
star matter. The deviation from the LS EOS case with 
K = 220 MeV, visible in Fig. 3, thereby arises only from 
the presence of pious and muons. 



Name 


K [MeV] 


Mmax/MQ 


Van [MoV] 


VsN [McV] 


VsN [McV] 


Vyy [MoV] 


180BG 


180 


1.15 


-28.2 


16.8 


-24.3 


0.0 


220BG 


220 


1.70 


-26.6 


28.5 


-22.8 


-38.0 


220g2.8 


220 


1.93 


-29.6 


65.7 


-23.0 


- 55.1 


220g3 


220 


1.95 


-26.8 


73.0 


-15.3 


- 10.3 


220pm 


220 


1.94 


-26.8 


24.1 


-24.5 


-10.3 



TABLE I: Properties of the different equations of state discussed in the text. The cold neutron star maximum mass is given 
for a non-rotating spherical star. 



Name 


ttAN 


Can 


ttBN 


Chjv 


asN 


cejv 


b-EN 


ttYY 


cyy 


6ee 


7 




MeV fm^* 


McV fm^*^ 


McV fm'"* 


McV fm^*^ 


McV fm^ 


McV fm^*^ 


McV fm^ 


McV fm^ 


McV fm 


MeV fm^ 




180BG 


-505.2 


605 


-434.4 


520.1 


10 


175 


214.2 











4/3 


220BG 


-340. 


1087.5 


-291.5 


9.32.5 


1.30. 


300 


214.2 


-486.2 


1553.6 


430 


2 


22()g2.8 


-270 


2:i()() 


-170 


21)00 


."lOO 


200 


211.2 


-100 


l.lOO 


i:!0 


2.8 


220g3 


-270 


4000 


-170 


2900 


450 


250 


214.2 


-90 


1000 


430 


3 


220pm 


-270 


4000 


-240 


3400 


130 


800 


214.2 


-90 


1000 


430 


3 



TABLE IL Parameter values of the different equations of state discussed in the text. 



V. COLLAPSE OF A COLD NEUTRON STAR 
TO A BLACK HOLE 

The goal of this Section is to show the numerical us- 
ability of the EOS, by implementing cold EOSs presented 
in Sec. IV B into a numerical code and studying the col- 
lapse of a cold neutron star to a black hole. It is not fully 
relevant to try and measure the effects of the presence of 
hyperons or additional particles in this case because it is 
difficult to find comparable physical settings: from Fig. 3 
it is clear that the maximal masses with or without hy- 
perons are quite different. Therefore, initial (unstable) 
neutron star models used as initial conditions are too 
different with different EOSs to provide comparable col- 
lapses to a black hole. The comparison can be performed 
in the case of a stellar core-collapse, starting from the 
same main-sequence massive star initial data and follow- 
ing then the collapse, bounce, stalling of the shock and 
collapse of the proto-neutron star to a black hole, with 
different EOSs at finite temperature. As this is beyond 
the scope of the current paper, it shall be the subject of 
a forthcoming study. 

A. Transition to a quark phase 

The physical model that we study here; is the collapse 
to a black hole of an unstable neutron star, i.e. with 
a central density higher than that corresponding to the 
maximal mass. During the collapse, the densities reached 
inside the neutron star can be much higher than several 
times the saturation density (see Fig. 4 and [68]). As the 
LS model for the EOS is non-relativistic, it can in prin- 
ciple allow for a sound speed greater than the speed of 



light. With the numerical code we are using (see Sec. VB 
hereafter), such a situation can lead to spurious oscilla- 
tions and instabilities destabilizing the whole simulation. 
This happens in the LS EOS or in most of the extended 
EOSs for high densities, typically beyond 5-10 fm~^, that 
is at more than twenty times nuclear matter saturation 
density. At these high densities, the nature of matter 
is anyway far from being well known. It seems rather 
natural to assume that there is a transition to a quark 
matter phase at some density and this is what we do here. 
Therefore we have used in our EOS a transition to a sim- 
ple model for quark matter, the MIT bag model [69] with 
massless quarks, implying a supraluminal sound speed. 
The transition is constructed using a Maxwell construc- 
tion. Within this very simple model the density for the 
transition can be adjusted by tuning the value of the bag 
constant. Since the aim of our paper is not to investigate 
the transition to quark matter (sec e.g. [16] for such a 
study in the context of core collapse), we try to push the 
transition density as far as possible without having super- 
luminal sound speeds. This means that we choose values 
of the bag constant much above commonly used values, 
leading to appearance of quark matter at a density of 
about ns = 2 fm~'^, the exact density thereby depends 
on the specific parametrization of the EOS used. 

With the quark phase transition presented here for the 
zero-temperature EOS, the sound speed always remains 
lower than the speed of light, suppressing all the possible 
instabilities coming from superluminal characteristics. In 
principle this problem could have been avoided by other 
prescriptions. For example, a causal EOS, i.e. a maxi- 
mally stiff EOS has been used in this context, too. We, 
however, think that a transition to quark matter better 
reflects the physics, although it should be kept in mind 
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that in this high density region we do not know much 
about the real physics and that our quark model is a 
very simplistic one. 



B. Model setup 

We use a gcncral-rclativistic hydrodynamics code Co- 
CoNuT [70], solving the continuity equation and the 
relativistic Euler equations for a perfect fluid by means 
of high-resolution shock-capturing methods [71]. Ein- 
stein equations, in isotropic gauge and maximal slicing 
are solved on a different grid with spectral methods [72]. 
This code, although able to perform 3-dimensional sim- 
ulations is run in spherical symmetry since we are only 
interested here in illustrating that the different EOSs can 
be successfully used in simulations of a collapse to a black 
hole. The EOS is tabulated, read from a file and in- 
terpolated for every grid point. Static, spherically sym- 
metric, neutron star initial data are obtained with the 
same EOS as the one used for the evolution, and with 
the same gauge, using the non-rotating version of the 
code described in [73] and the LORENE library [74]. The 
initial star is chosen to be a spherical configuration on the 
unstable branch, with a gravitational mass decreasing if 
one increases central density. An unstable neutron star 
is not bound to collapse to a black hole: it can expand 
in order to "migrate" to the stable branch, i.e reach a 
lower density with the same baryon mass, such that it is 
stable with respect to radial oscillations (see also [75]). 

The star is then perturbed by amplifying the radial 
density profile by one percent. This procedure ensures 
that the star will collapse to a black hole and not mi- 
grate to the stable branch. The standard picture of such 
a numerical model is that the star collapses until general- 
relativistic effects become dominant (see e.g. [68] for a 
complete description). Among these is the "frozen-star" 
effect which comes from our choice of time gauge (max- 
imal slicing), avoiding the appearance of a central sin- 
gularity. It implies that several quantities do no longer 
evolve near the center of the star, whereas some of the 
metric coefficients develop huge gradients, limiting the 
simulation in time. One sign that evolution is frozen near 
the center of the star is given by the fact that one metric 
coefficient, the lapse a, representing the ratio between 
the physical time measured by the Eulerian observer and 
the coordinate time, is tending toward zero. The second 
general relativistic effect is the formation of an appar- 
ent horizon at a finite distance, growing further outward 
until including all the neutron star matter. This is an 
evidence of the formation of a black hole. The 2-surface 
defining the horizon is tracked in our code by an appar- 
ent horizon finder [76] , which enables us to compute the 
baryon mass inside the black hole, too. 
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FIG. 4: Profiles of central density (solid line) and central 
lapse (dashed line) as functions of time during the collapse 
of a perturbed unstable neutron star to a black hole. Top 
figure was obtained with the LS EOS [3] with K = 220 MeV, 
the two others by EOSs described in Tabs. I and II. Vertical 
dotted lines give the time of the formation of the apparent 
horizon. 



C. Results of simulations 

We have run our code on six different tabulated EOSs 

displayed in the left column of Tab. III. In this table, 
LS180 and LS220 stand for the Lattimer-Swesty EOS 
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EOS 


tAH [rns] 


tgg [ms] 


ng- [fm--^] 


AA4b 


AMg 


LS180 


0.223 


0.23 


18.6 


2 X 10"*' 


4 X 10"-^ 


180BG 


0.188 


0.216 


114 


7 X 10"** 


2 X 10"^ 


LS220 


0.209 


0.216 


6.19 


2 X IQ-® 


4 X IQ-^ 


220g2.8 


0.387 


0.396 


20.8 


5 X 10"'^ 


6 X 10"^ 


220g3 


0.327 


0.334 


9.51 


6 X lO"'' 


5 X 10"^ 


220pm 


0.215 


0.22 


3.24 


IQ-^ 


6 X IQ-^ 



TABLE III: Characteristics of the six collapses to a black 
hole studied here. EOS names arc detailed in the text, Iah 
is the apparent horizon formation time (since the starting of 
the collapse), tgg the time at which 99% of the baryon mass 
has gone into the black hole, ng''" is the maximum central 
density reached during the collapse, AMb and AMg are the 
relative conservations of baryon and gravitational masses, re- 
spectively. 



dition of hyperons makes the EOS even softer and, at 
the end the collapse vi^ith EOS 180BG, it reaches such 
high densities that they no longer seem realistic (beyond 
100 fm~^). Finally, Tab. Ill also gives some error indica- 
tors: AMb is the relative variation of the baryon mass 
(number of particles times their rest mass), AMg is the 
variation of the system's gravitational mass (as deduced 
from the asymptotic behavior of the gravitational field) . 
The conservation of baryon mass is directly imposed by 
solving the equation for the conservation of the baryon 
current, whereas the conservation of gravitational mass 
is only an indirect consequence of the solved equations 
and of the spherical symmetry (no gravitational waves). 
Therefore, this last indicator is a good estimate of the 
overall accuracy of a run. From all these results, we can 
claim that, apart from EOS 220BG and 180BG, the cold 
EOSs d(;rivc!d here are suitable for numerical simulations 
in the demanding model of the collapse to a black hole. 



with the incompressibility K = 180 and K = 220 MeV, 
respectively. The properties of the four other EOSs are 
given in Tabs. I and II. The EOS labeled 220BG in these 
tables could not give any reliable result in the simula- 
tions, because too much numerical noise appeared al- 
ready in the initial data. Neutron star collapses with 
any of the six EOSs listed in Tab. Ill would lead to the 
formation of a black hole in a time of a few tenth of a 
millisecond. Time evolution profiles for the density and 
the central value of the lapse a are given in Fig. 4 for 
EOSs LS220, 220pm and 220g2.8. On each of these plots 
the central density is increasing, eventually with some 
oscillations coming from the focusing of the initial per- 
turbation, before reaching a maximum value and then 
decreasing. This decrease is mostly due to the finite res- 
olution at the star center, as the density should tend 
toward a given value (frozen-star picture). This spuri- 
ous decrease is not an issue because it appears after (or 
about the same moment as) the formation of the appar- 
ent horizon and therefore in a region inside the black 
hole that cannot influence the matter still falling onto it. 
This problem can in principle be cured using the exci- 
sion technique (removing a neighborhood of the center, 
replacing it by boundary conditions [77]), which we plan 
to implement in our code in the near future. 

On Fig. 4, the lapse is decreasing toward zero, as ex- 
pected, and the apparent horizon forms sufficiently soon 
to let most of matter remaining outside the black hole to 
enter it. The simulations are ended when the radial gra- 
dients of the gravitational fields become too high to be 
well described with spectral methods. As seen from the 
values of tgg Tab. Ill, this happens after all matter has 
been swallowed by the apparent horizon and therefore 
one is left with a static Schwarzschild black hole, with no 
evolution outside the horizon. From Fig. 4 and Tab. Ill, 
one can notice that the maximal density reached during 
the collapse strongly depends on the EOS. In particular, 
the EOSs with an incompressibility modulus K = 180 
MeV being softer, matter is more compressed. The ad- 



VI. RESULTS AT FINITE TEMPERATURE 

Let us now discuss the behavior of the EOS with the 

different parametrizations at finite temperature. As men- 
tioned earlier, simulations of core collapse with massive 
progenitors show that rather high temperatures of sev- 
eral tens of MeV and even more at proto-neutron star 
densities are reached, see e.g. [17, 78]. It is clear that the 
thermal energy is in favor of the production of additional 
particles such as hyperons and mesons or nuclear reso- 
nances. The abundances of hyperons, pions and kaons 
measured from heavy ion collisions indeed indicate that 
they are produced during the collisions. One should of 
course insist on the fact that the conditions in heavy 
ion collisions are different from core collapse events: the 
baryon densities are lower and the temperatures are prob- 
ably slightly higher, the timescales are such that no weak 
equilibrium for strangeness is achieved and the difference 
in neutron and proton densities is much less pronounced 
than in core collapse events or in neutron stars. But the 
results show clearly that the thermal production of those 
additional particles is important. 

At finite temperature and without assuming /3- 
equilibrium, the EOS is a function of three variables 
which are generally chosen to be T, ns, Vg- We now dis- 
cuss the properties of the EOS as a function of these vari- 
ables. Of course we cannot cover the whole range, but 
we shall choose some particular conditions. In Fig. 5 we 
show the different hyperon fractions as functions of elec- 
tron fraction for a temperature of T = 60 MeV and two 
different densities, ub = 0.15 fm''' and ub = 0.3 fm~^. 
In the simulations of Sumiyoshi et al. [17] the first den- 
sity corresponds to the conditions of a collapse with a 
4OAf0 progenitor at bounce at about 10 km from the 
center. Of course the exact thermodynamic conditions 
in a simulation depend on the EOS, so that this is just 
to say this is a typical situation within a proto-neutron 
star after bounce. High densities are typically reached 
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FIG. 5: (Color online) Fractions of different hyperons as a function of the electron fraction at a temperature of 60 MeV and for 
ub ~ 0.15 fm~'^ (right) and n_g = 0.3 fm~'^ (left) corresponding roughly to once and twice nuclear matter saturation density. 
The fractions of S"'^, are not shown since they are below one percent. 



in the post-bounce phase. We only show the fractions of 
A, and S^-hyperons since the fractions of the other 
hyperons are very small, between one permille and one 
percent. 

At hb ~ 0.15 fm^"^, the A and fractions show 
clearly two groups: the parametrizations without pions 
and muons and the one including them. The reason for 
the different behavior of the EOS if pions and muons are 
included is mainly due to the presence of muons. This 
can be understood rather easily. Under the present con- 
ditions muons are almost equally abundant as electrons, 
so that in order to conserve global charge neutrality, the 
hadronic charge fraction for a given is increased by 
almost a factor of two with respect to matter without 
muons. This is clearly reflected in the curves. Let us 
stress, however, that the muon fraction is determined 
by the degeneracy factor /i;^/T so that it depends on 
the ratio of chemical potential and temperature. The 
curves are calculated assuming the same chemical poten- 
tial for electrons and muons. In the case that neutrinos 
can freely leave the system, we have fi^- — ^(.- — l^q: 
but in the hot proto-neutron star neutrinos are trapped 
so that there are nonzero lepton number chemical poten- 
tials which are not necessarily the same for electrons and 
muons as mentioned already earlier. In general electron 
neutrinos are the most abundant ones, so that assum- 
ing an equilibrium Fermi-Dirac distribution the lepton 
number chemical potential for electrons is expected to 
be higher than that of the other leptons. We therefore 
expect that in a realistic simulation the chemical poten- 
tial for muons is lower than that for electrons reducing 
the number of muons. It is, however, difhcult to esti- 



mate quantitatively the reduction. We can only say that 
we think that our results represent an upper limit for the 
importance of the effects muons can have on the EOS. 

The E~ fraction shows a stronger dependence on the 
hyperonic interaction. Remind that the two parametriza- 
tions 220g3 and 220g2.8 contain a rather strong repulsion 
for the single-particle potential, see Tab. I, whereas 
the repulsion is much less strong for the others. This 
observation explains the reduction of the -fraction in 
these two models with respect to the others. The overall 
increase in the -fraction with decreasing is the well 
known effect that in neutron rich nuclear matter nega- 
tively charged hyperons are favored. To a smaller extent 
the same is observed for the S~, which is, however, less 
abundant overall due to its higher mass. Altogether the 
hyperon fractions reach maximum values of about 2-8 %. 

In Fig. 6 the pressure, the sound speed and the internal 
energy per baryon with respect to the proton mass rrip are 
displayed for the same densities and temperature. The 
latter quantity is defined as 

e = ^— - 1 , (9) 
riBrrip 

with e denoting the total energy density. These three 
quantities are key ingredients for the hydrodynamic 
simulations. We show the results for the different 
parametrizations of the hyperonic interactions as well as, 
for comparison, the LS EOS with K = 180 McV and K 
= 220 MeV, too. 

We again observe that the EOS including muons be- 
haves differently from all the others. The reason has been 
explained above. In particular, the usual softening of the 
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FIG. 6: (Color online) Thermodynamic quantities as a function of the electron fraction at a temperature of 60 MeV and for 
ns = 0.15 fm~'^ (right) and n_g = 0.3 fm~"^ (left) corresponding roughly to once and twice nuclear matter saturation density. 
The upper panels show the pressure, the middle ones the internal energy per baryon with respect to the proton mass (see text 
for the definition) and the lower ones the sound speed squared. 



EOS from additional degrees of freedom by adding dif- 
ferent types of particles is not seen. The point is that 
this softening is overcompensated by the effect of the in- 
creased hadronic charge fraction induced by the presence 
of muons, see above. For the other EOS including only 
hyperons the softening is indeed seen, as expected it is 
more pronounced at higher density (left panels). Com- 
pared with the purely nuclear EOS e is higher including 
hyperons. This is not very surprising neither since the 
hyperons are more massive than nucleons and therefore 
replacing a nucleonic state with a hyperonic one in gen- 
eral increases the energy density. The modifications of 
pressure and energy density due to the presence of addi- 
tional particles is reflected in the sound speeds, too. As 
already mentioned above, the hyperon fractions increase 
with decreasing electron fraction, so that the effects on 
the thermodynamic quantities increase with decreasing 
Ye, too. 

At UB = 0.15 fm~^ (right panels) the interaction 
has only little effect on the thermodynamic quantities, 
three groups of EOS can clearly be distinguished: the 
two LS ones, those with hyperons and the EOS includ- 
ing in addition pious and muons. The differences aris- 
ing from the different interaction between the two LS 
ones, and between the four hyperonic ones are only very 
small. This can be understood since at high tempera- 
ture and low density the kinetic energy should dominate 
and the interaction terms should in turn be less impor- 



tant. At higher densities, as can be seen from the figures 
at ub = 0.3 fm~^ (left panels), indeed the infiuence of 
the interaction is stronger. The difference in pressure be- 
tween the two LS ones is of the order of 10% and between 
the four hyperonic ones of the order of 15%. It should 
be stressed, however, that the softest one, 180BG, gives 
a maximum mass for a cold spherical neutron star be- 
low 1.2 Mq and that it is thus not very realistic, see 
Sec. IV B. The difference between the three remaining 
hyperonic EOS is much smaller. 

At which temperatures the additional particles in the 
EOS start to play a role? In order to answer this question 
we display in Fig. 7 the fractions of A, , and S~ and in 
Fig. 8 the pressure, e and the sound speed as a function 
of temperature. The densities are the same as before, 
ub — 0.15 fm""^ on the right and hb — 0.3 fm^'^ on the 
left. We have chosen a relatively low electron fraction. 
Ye — 0.1, because we want to show an upper limit case, 
i.e. an optimistic estimation of the effect of the additional 
particles on the EOS. At the smaller density hyperons ap- 
pear at about 25 MeV, independently of the EOS used. 
The first one to appear is the A-hyperon. As expected the 
hyperon fractions rise with temperature, reaching about 
10% for A and and 5% for E!~ at a temperature of 100 
MeV. The differences in the parametrization of the hy- 
peronic interaction does not induce large differences. Let 
us make, however, two remarks. First, the particular fea- 
tures of the EOS with pions and muons have already been 
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explained and it is thus clear why the hyperon fractions 
are systematically lower for this EOS than for the others. 
Second, for the S^-fraction, the two paramctrizations 
by Balberg and Gal, 180BG and 220BG, clearly show a 
higher E~ -fraction than the other ones. The reason is 
again the strong E~ repulsion in the parametrizations 
220g2.8 and 220g3. 

At ns = 0.3 fm^'^, hyperons appear at lower tem- 
peratures, depending on the EOS between roughly 15-25 
MeV. The differences between the EOS are more pro- 
nounced than at lower density. As mentioned earlier, 
this can be imderstood from the fact that, for sufficiently 
low density, due to thermal effects the kinetic energy 
should be dominant rendering the details of the inter- 
action less important. For this higher density the A- 
hyperon is clearly the most abundant one, attaining be- 
tween 10 and 17 % at T = 100 MeV. The S^-fraction at 
this temperature lies between 4 and 8% and for the E~, 
slightly less abundant, the fraction reaches 2-6 %. 

The thermodynamic quantities, in particular pressure 
and sound speed clearly show the appearance of hyper- 
ons, inducing a softening in the EOS. For the lower den- 
sity, the modifications in the pressure due to the pres- 
ence of hyperons stays, however, relatively small up to 
T = 100 MeV, whereas e and the sound speed show more 
important deviations between the purely nuclear case and 
the different cases with additional particles. Again the 
particle content of the EOS has more infiuence on the 
behavior of the thermodynamic quantities than the de- 
tails of the interaction. This is, as already noticed be- 
fore, not true at twice this density, where the different 
parametrizations give different results for pressure and 
sound speed. For the energy density the differences are 
smaller. 

We recover most of the features discussed up to now as 
a function of density, too. This can be seen from Fig. 9, 
where the hyperon fractions are shown and from Fig. 10, 
where the thermodynamic quantities are displayed for 
Fe = 0.1. On the left panels T = 25 MeV, on the right 
panels T — 40 MeV. An interesting point which we have 
not seen before because the density has been too low is 
that at about 2.5no a transition takes place, strongly in- 
creasing the hyperon fraction and with a strong effect on 
the thermodynamics. The A-fraction, for instance, can 
in certain models be larger than 30 %. The thermody- 
namic quantities, in particular the pressure, reflect this 
transition. A thorough discussion of the thermodynamics 
of this transition togethcir with a detailed analysis of the 
stability is of order but beyond the scope of the present 
paper. 

From the EOS alone we, of course, cannot answer 
the question whether the modifications in the thermody- 
namic quantities due to the presence of hyperons, pions 
and muons are relevant for the dynamics of a core col- 
lapse event or a neutron star merger. Here, we do not 
want to try to answer this question, but in order to get an 
idea we compare the pressure, energy density and sound 
speed profiles for a hot proto-neutron star. The data for 



this profile, thus T,nB,Ye, are shown on the right panels 
of Fig. 11 as a function of the radius. They are issued 
by a ID simulation of the collapse of a 40 Mq progen- 
itor with full Boltzmann neutrino transport employing 
the LS EOS with K = 180 MeV, see Ref. [79, 80], at 
about 400 ms after bounce. The left panels show the 
pressure, e and the sound speed as a function of the radius 
for the different EOS corresponding to the given values 
of temperature, baryon density and electron fraction at 
this radius. Of course, this procedure does not give cor- 
rect proto-neutron star profiles since these depend on the 
EOS. However, we think that with this remark of caution 
in mind, the comparison of the profiles is interesting and 
can give hints on the importance of the modifications in 
the high density and high temperature part of the EOS. 

Let us first examine the data. The temperature is 
about 40 MeV at the center, rising to more than 80 
MeV at about 10 km from the center and decreasing 
then rapidly to a value between 5-10 MeV. The density 
is maximal at the center with a value slightly above 3no 
decreasing to below saturation density at about 10 km 
from the center. has a value of about 0.3 at the cen- 
ter. It decreases until about 15 km, where the value is 
about 0.1 and rises then, reaching 0.5 at about 30 km 
from the center. 

From the behavior of T^ns and Yg we would expect 
sensible modifications of the thermodynamic quantities 
only within a radius of about 15 km from the center of 
the proto-neutron star. This is indeed the case as can be 
seen from the left panels. Let us start with comparing 
the two different versions of the LS EOS. In the region 
very close to the center, up to roughly 8 km, the pres- 
sure and the sound speed show differences between the 
two LS EOS. This is understandable since only in this 
region the density is high enough to allow for the differ- 
ences in the nuclear interaction to play a significant role. 
Below saturation density at low temperatures the nuclear 
EOS is relatively well constrained so that in this region 
it would be surprising to see large differences. For high 
temperature and low density the kinetic part becomes 
dominant, so that no large differences due to the details 
of the nuclear interaction are to be expected. Therefore 
only the high density region remains where different nu- 
clear EOSs can show very different behavior. Remember 
that we are not interested here in the details of the nu- 
clear composition, which can have an influence on the 
thermodynamics and on the dynamics of a core collapse 
event too, see [5]. These occur mainly below saturation 
density and temperatures below 10 MeV, so that they 
cannot be resolved on the scales we are examining here. 

Concerning the comparison with the extended EOS, we 
can see differences up to a radius of about 15 km. These 
differences are in general more pronounced than those 
between the two version of the LS EOS. The pressure 
at the center varies between 60 and 100 MeV/fm'^ from 
the "softest" to the "stiffest" EOS. The lowest pressure 
is obtained for 180BG, the second lowest for the LS EOS 
with K = 180 MeV and all the others give values above 
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FIG. 9: (Color online) Same as Fig. 5 but as a function of baryon number density for an electron fraction Fe = 0.1 at a 
temperature of 25 MeV (left) and 40 MeV (right). 




FIG. 10: (Color online) Same as Fig. 6 but as a function of baryon number density for an electron fraction Ye ~ 0.1 at a 
temperature of 25 MeV (left) and 40 MeV (right). 
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FIG. 11: (Color online) Pressure, energy density and sound speed (left) for different EOS as a function of radius for a proto- 
neutron star profile (right) about 400 ms after bounce with a 40 Mq progenitor. Data (right) are from a simulation by IVI. 
Liebendorfer employing the LS EOS with K = 180 MeV. 



90 MeV/inr^. This means that the influence of the addi- 
tional particles on the pressure in this proto-neutron star 
is less important than the value of the nuclear incom- 
pressibility. We of course expect this conclusion to no 
longer hold if the central density, only about 3no at the 
stage we are examining here, increases. And of course, as 
mentioned before, we have to be careful since we should 
recalculate the proto-neutron star profile with having a 
new EOS. 

For e we clearly have three groups of EOSs, the LS 
ones, those with hyperons and that with hyperons, pi- 
ous and muons. The differences are most pronounced at 
the temperature maximum. This latter point is related 
to the fact that the energy density strongly depends on 
temperature. Moreover, the profile in e closely follows 
the temperature profile. The previous discussion has al- 
ready shown why e has a distinct behavior depending on 
the particle content. The sound speed reflects the differ- 
ences in pressure and energy density and can vary at the 
center between roughly 0.45c and 0.6c. 

The results for 220pm shown in Fig. 11 have again 
been calculated assuming the same chemical potential 
for electrons and muons. As discussed above, this rep- 
resents probably an upper limit for the muonic effects 
on the EOS. Here, we have in principle the muon neu- 
trino fraction at hand. Thus, assuming that the muon 
neutrinos are in thermal equilibrium, i.e. that they are 
described by a Fermi-Dirac distribution, we can deter- 
mine the muon lepton chemical potential, and from this 



the muon chemical potential. We have computed the pro- 
files with this muon chemical potential and, as expected, 
the differences are smaller, but the general trends are the 
same and in particular the results for 220pm are still very 
different from the other EOSs. 



Up to now we have assumed that weak equilibrium 
with respect to strangeness is achieved. If, on the con- 
trary, we assume that we have no weak strangeness 
changing reactions, which would correspond to having 
reaction timescales much longer than the hydrodynamic 
timescale of 10~^ s, strangeness becomes a conserved 
quantum number. We do not consider this as a realis- 
tic scenario since the timescales estimated for the rel- 
evant processes are of the order of 10~^ or below, see 
e.g. [81]. We, however, find it instructive to compare 
our results with this extreme case. The hyperons have 
all negative strangeness, so that populating hyperonic 
states leads to a net negative strangeness. Typical pro- 
duction reactions for hyperons via the strong interaction, 
for example n + n — ^ n + K + A, are strangeness conserv- 
ing and kaons are produced with positive strangeness. 
Thus kaons are the natural candidates for assuring van- 
ishing net strangeness in thermal equilibrium. However, 
their mass of niK ~ 500 IVIeV is rather high compared 
with the relevant chemical potentials and temperatures 
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so that they are not very abundant ^ if they are consid- 
ered as an ideal gas. This in turn strongly suppresses the 
hyperon fractions and therefore the effects on the EOS 
compared with the more realistic scenario of strangeness 
changing weak equilibrium. 



VII. SUMMARY 

At densities above roughly nuclear matter saturation 
and temperatures above several tens of MeV, an equa- 
tion of state based uniquely on nucleonic degrees of free- 
dom and electrons is no longer realistic since many other 
states will appear. We have presented here an extended 
version of the LS EOS [3] including as additional parti- 
cles hypcrons, pions and muons intended to improve on 
the high density and high temperature part. For zero- 
temperature high density matter this question has al- 
ready been studied for many years but up to now only 
very few work exists for finite temperature. The main 
problem in this type of exercise is that the interaction, 
which is already not well known for nucleons, is even less 
known for hyperons. We have adapted here a very sim- 
ple phenomenological approach based on the hyperonic 
model of Balberg and Gal [46]. The parameters of the 
model have been readjusted in order to be compatible 
with available hyperonic data and in particular the ob- 
servation of an almost two solar mass neutron star, PSR 
1614-2230 [22]. Taking these constraints into account 
there still remains some freedom, so that we have dis- 
cussed several parametrizations of the EOS in order to 
get an idea of the uncertainty. The ultimate goal should 
be, of course, to have a reliable microscopic approach 
to hyperonic matter compatible with data, but awaiting 
this step, we can, phenomenologically, study the effect of 
these additional particles on the thermodynamics of the 
system. The results show that key thermodynamic quan- 
tities as pressure, energy density and sound speed are 
influenced by the additional degrees of freedom in a non- 
negligible way. The threshold temperature for the ap- 
pearance of hyperons at saturation density lies at about 
25 MeV, depending on the particular model applied. Due 
to the shift in the hadronic charge for a given electron 
fraction, muons seem to strongly influence the EOS in 
the regions where they become abundant. 

We have concentrated here mainly on hyperons as ad- 
ditional particles. Nuclear resonances have not been con- 
sidered for the moment. A study of this point is kept 
for future work. In another respect our EOS could be 
improved: we should not treat pions as a free gas but 
include interactions. At high temperatures this should 
not be very important, but for the low temperature and 



^ Note that we are not discussing here possible medium modifi- 
cations of the kaon properties which could lead to higher kaon 
abundances. 



high density regime we expect it to have some influence. 

The main application of our EOS should be astrophys- 
ical systems, first of all core collapse events of massive 
progenitors, collapsing eventually to a black hole. Neu- 
tron star mergers could be another application. We have 
demonstrated that our EOS can be successfully used in 
a numerical simulation of the collapse of a cold neutron 
star to a black hole. Of course the effect on realistic sim- 
ulations including finite temperature has to be tested. 
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Appendix A: Some technical issues on the 
construction of the EOS tables 

Complete EOS tables, including the extended versions, 
will be prepared and made publicly available. In this 
Appendix we would like to discuss some technical issues 
encountered in the construction of those tables from the 
LS original routine. We remind that we have modified 
the original routine in two respects: (1) correcting the 
binding energy of the a-particle (see Fig. 1); (2) extend- 
ing the routine at densities ^ 10~^ fm~^, by extending 
the validity of the Maxwell and boundary construction 
files. This also permits to veriiy a correct matching with 
the low density EOS. However, some convergence issues 
at low temperatures and proton fractions or near criti- 
cal temperature and density still remain. Convergence 
problems manifest in two ways: (1) no solution of the 
equilibrium equations, Eqs. (3.2) in Rcf. [3], is found (es- 
pecially at low temperature and electron fraction [31]), 
(2) the solution is discontinuous with respect to adjacent 
points in density, temperature, and/or electron fraction. 
In the latter case, the discontinuities appear as: (i) a 
rapid changes of the regime (with/ without nuclei), which 
comes from a cross over the boundaries, (ii) a conver- 
gence towards a point far from the adjacent ones. The 
first kind of pathology might be due to the fact that a 
phase coexistence (and not a phase mixing) is considered 
to model the phase transition. The second one might be 
due to the sensitivity of the solution with respect to the 
starting point values in the minimization routine. 

The construction of a table partially overcomes these 
discontinuities, since, when looking up the table, an in- 
terpolation is done and "critical" points in-between the 
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grid points are usually avoided. 

Another point to bear in mind when constructing the 
EOS tabic is the speed of sound. This is a crucial quan- 
tity in the hydrodynamic simulations, since it regulates 
the speed at which the information is propagated and 
it determines the time step in finite difference explicit 
schemes. Throughout the calculation, due to convergence 
problems, it might appear that the speed of sound is ei- 
ther less than zero or superluminal. This is of course non 
physical, so for those (rare) points we have recalculated 
the EOS replacing the value of the speed of sound with 
the one obtained in the case of a Fermi gas (in units of 
the speed of light) (see e.g. [82]): 



{hey 



3(mec2)^ 



(Al) 



This replacement does of course not concern the region 
at high densities where the sound speed becomes super- 
luminal due to the non-relativistic character of the LS 
EOS, see Section VA. 



We have compared our tables to the O'Connor and 
Ott ones [83]. However, even if there is an agreement 
in the range where the LS routine is employed, it is not 
straightforward to make an exact comparison since the 
nuclear parameters that they employ are slightly different 
from the ones used in this paper. 
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